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ABSTRACT 

Context. High-contrast imaging is currently the only available technique for the study of the thermodynamical and compositional 
^ ' properties of exoplanets in long-period orbits, comparable to the range from Venus to Jupiter. The SPICES (Spectro-Polarimetric 

CIh, Imaging and Characterization of Exoplanetary Systems) project is a coronagraphic space telescope dedicated to the spectro- 

polarimetric analysis of gaseous and icy giant planets as well as super-Earths at visible wavelengths. So far, studies for high-contrast 
imaging instruments have mainly focused on technical feasibility because of the challenging planet/star flux ratio of lO^^^-lO"'" re- 
quired at short separations (200 mas or so) to image cold exoplanets. However, the main interest of such instruments, namely the 
^ I analysis of planet atmospheric/surface properties, has remained largely unexplored. 

Aims. The aim of this paper is to determine which planetary properties SPICES or an equivalent direct imaging mission can measure, 
considering realistic reflected planet spectra and instrument limitation. 

Methods. We use numerical simulations of the SPICES instrument concept and theoretical planet spectra to carry out this performance 
^ \ study. We also define a criterion on the signal-to-noise ratio of the measured spectrum to determine under which conditions SPICES 

\ j^ , can retrieve planetary physical properties. 

fsj ' Results. We find that the characterization of the main planetary properties (identification of molecules, effect of metallicity, presence 

QQ ' of clouds and type of surfaces) would require a median signal-to-noise ratio of at least 30. In the case of a solar-type star <10 pc, 

^s^ I SPICES will be able to study Jupiters and Neptunes up to ~5 and ~2 AU respectively, because of the drastic flux decrease with 

• , separation. It would also analyze cloud and surface coverage of super-Earths of radius 2.5 Earth radii at 1 AU. Finally, we determine 

C^ ■ the potential targets in terms of planet separation, radius and distance for several stellar types. For a Sun analog, we show that SPICES 

^^ ' could characterize Jupiters (M > 30 Earth masses) as small as 0.5 Jupiter radii at <2 AU up to 10 pc, and super-Earths at 1-2 AU for 

Cn I the handful of stars that exist within 4-5 pc. Potentially, SPICES could perform analysis of a hypothetical Earth-size planet around 

a Cen A and B. However, these results depend on the planetary spectra we use, which are derived for a few planet parameters assum- 
ing a solar-type host star. Grids of model spectra are needed for a further performance analysis. Our results obtained for SPICES are 
also applicable to other small (1-2 m) coronagraphic space telescopes. 
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1. Introduction from radial velocity surv eys ([Mayor etaP 201 1[ hereafter RV) 

and transit su rveys (e.g.. | Leg er et al.' 2009; Ch arbonneau et all 

_, , £,.•,.-• , • u jj- 2009; Batalha etaflboUriBorucki et al. 2012). There are so 

The exoplanet held in astrophysics IS extremely rich and diverse. ' • i j j j i i j 

„ j^.-.u . ■ .■ . t ■ T. ■ many planet categories already detected even with the detec 

From detection to characterization, many techniques are being ■ , ■ r , , a ,r^xr ■ ■ 

used or developed to address the fundamental questions about 

planetary formation and evolution. Exoplanets span a number 

of categories much larger than the Solar System's planets do. 

Since the first discovery, several unexpected types of exoplanets „ ■ ,, • ■^- -^ . i ■ j 

„ ,„ ,u u . t •. .-KK o /-v 1 inrvc, u- u Current eliorts are to expand their sensitivity to longer-peHod 

werefoundfromthehot Jupiters (May or & Oueloz 1995) which . „ , , atr:, — „ Z , — i ^rvrv-rl fe o -r^ — ^n 

.r. .,.r., .I.e. t. th.ir hr^.^.^.r. ( ^(^ H^ ATH t. th. ...ni.ti.n ^nd Smaller cxoplancts ( iUdry & Santos|2007HSeager & DemingI 



tion biases of the current methods (RV, transits, imaging, mi- 
crolensing) that several instruments/missions will be needed to 
cover the whole field. Methods such as RV and transits ap- 
pear to be effective at probing for large close-in exoplanets, and 



are very close to their host stars (<0.05 AU) to the population 
of super-Earths (massive telluric planets) which starts to emerge 



2010) ■ The detection and characterization of long-period/wide- 
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orbit (>1 AU) planets are, however, still difficult to accomplish. 
On the one hand, it requires long-duration monitoring with very 
stable instruments to detect them from a RV or transit survey. 
On the other hand, direct imaging has to tackle the high contrast 
at small angular separation that exist between the planet and its 
host star. 

An extrapolation of the period distribution of giant planets 
discovered by RV surveys suggests that a large population of 
these objects at separations la rger than 5 AU still remains to 
be revealed ( Marcv et al.ll2005h . These planets, in particular the 
population between 5 and 20 AU, are very important for con- 
straining theoretical mechanisms of planetary formation, since 
it may reside at the bo undary between the core-accret ion and 
disk instabihty regimes (lAlibert et al .1120 lit lBossll201lh . When 
it comes to the atmospheric characterization of these planets, 
direct imagi ng will pro bably be the most productive te chnique. 
Since 2005 (IChauvin et al. 2005; Neu hauser et al.ll2005h . several 
massive giant planet candidates were imaged around young stars 
(<200 M yrs), the most emblematic being the fou r planets around 
HR 8799 (Marois et al. 2008:'Ma rois et 31.12010 *) and/? Pictoris b 
(JLagrange et al. 2 009, 2010) . We note that the pl anetary nature 
of the Fomalhaut companion ( Kalas et al] l2008l) has been re- 
cently questioned (Jan son et al.l2012h . Following these discover- 
ies, spectra we re obtained f or a few planets (e.g. , Mohanty et al. 
2007; iJanson et al.1 |2010; Patience eLsD l20ia IBowler et al.1 



20101). A first generation of instruments precisely optimized for 
the detection and spectral characterization at near- and mid- 
infrared (IR) wavelengths of young giant planets will see first 
light in the present dec ade: SPHERE (Beuzitetal. 2008), GPI 
dMacintosh et al.' '2008), HiCIAO (Hodappetal. 2008), P1640 
Phase II (.Hinklev et al...2011.) and FLAO ( Esposito et al,.2O10.) 
on ground-based telescopes, an d JWST ( Clampin] l2010h and 
SPICA Coronagraph Instrument (lEnvaet al.ii2011i) in space. In 
the next decade, planet finders on E xtremely L arge Telescopes, 
such as EPICS (lKasperetalJl2010h and PFI dMacintosh etTaP 
120061) . may oifer the ability to observe mature gas giants, ice 
giants and possibly super-Earths in the near-IR. Detailed stud- 
ies were performed to consider the feasibility of large aperture 
coronagraphs and large baseline interferometers for the detec- 
tion of Earth twins from space. These studies identified areas 
of technological development that need to be first addressed, 
which will delay the realization of Terrestrial Planet Finder mis- 
sions until ~2025-2030. Meanwhile, observations have demon- 
strated the extreme diversity of planets. This led to a growing 
consensus within the community that we will need to study 
all planet types in order to have a complete u nderstanding of 
their formation and evolution (ISchneider et al.l2 008). To address 
some parts of these questions, a family of small space missions 
(mainly coronagrap hs) have been proposed (Guvon et al. 2010b; 
iTrauger et al.ll2010l) for analyses of ice giants and super-Earths. 
The study of these missions has mainly focused on technical 
feasibility, as it i s a challenge to achiev e large contrast close 
to a bright star (iTrauger & Traubl l2007[ iGuvonetal.l l2010ai 
iBelikov et al.l l2010l) . Another area of study is the estimated 
number of observable exo planets of a given type (ITrauger et al.l 
|20 1 0*; 'Guvon et al.'20 1 ObV based on assumptions of their densitv 
distribution. Cahov et al. (2009) consider signal-to-noise ratios 
(SNRs) of broad-band (R = 5) photometric measurements as- 
suming planets with grey albedos. The main interest of such mis- 
sions, the atmospheric c haracterization, ho wever, has remained 
highly unexplored so far. ICahov et aP (12010 ) analyze colors and 
coarse spectra (R = 5 and R = 15) of Jupiter and Neptune atmo- 
sphere models, but without including instrument limitations such 
as throughput and noise. From the colors of Solar System plan- 



ets, iTraubl (l2003h suggested that for planets too faint for spec- 
trometry, even coarse colors could help to distinguish between 
planet types. Cahov et al. (2010) showed that while a color cri- 
terion could not be a means to uniquely distinguish between 
planet types, due to the intertwined contributions of factors such 
as metallicity and planet-star separation, colors would still pro- 
vide some constraints on possible planet types. 

Among all of these space coronagraph concept missions, 
the most recent, SPICES (Spectro-Polarimetric Imaging and 
Characterization of Exoplanetary Systems) was submitted to the 
ESA Cosmic Vision call for medium-class missions in 2010 by a 
consortium of Eur opean institutes with American and Japanese 
participation^ (B occaletti et alJl2012l) . SPICES has a twofold 
motivation: 1/the systematic atmospheric characterization of gas 
and ice giants as well as super-Earths in the solar neighborhood, 
and 2/ the development and validation of key technologies in or- 
der to prepare future direct imaging projects dealing with Earth 
twins spectral characterization. The science objecti ves and the 
technical concept of SPICES are described in Boc caletti et al.l 
(2012). The main science driver is the study of planetary systems 
as a whole for the understanding of planet formation and evolu- 
tion. With a maximum imaged field of view of ~13", SPICES 
will focus on targets previously identified using other methods 
(planets and circumstellar disks), but can also detect new plan- 
ets such as outer planets in known planetary systems and exo- 
zodiacal disks <100 zodis. A preliminary estimation of the num- 
ber of characterizable planets gives an order of magnitude of 100 
objects for an allocated time of three years over the five years 
of the mission (Boccaletti et al. 2012). The main purpose of the 
instrument is to obtain flux and polarization spectra at visible 
wavelengths of cold/mature exoplanets, especiall y those previ- 
ously discovered by RV surveys (lUd rv & San tos! l2007b or as- 
trometry with GAIA (ICasertano et al.l l2008). These surveys will 
provide the orbital elements and minimum mass of planets, but 
not the radius which determines the amount of reflected light 
together with the albedo. It is therefore essential to perform ac- 
curate measurements of spectra to possibly distinguish between 
planet types. 

In this work, we have developed a numerical simulation to 
model the instrument concept of SPICES. Under realistic as- 
sumptions of noise, instrument performance and reflected exo- 
planetary spectra, we test the ability to distinguish between spec- 
tra of planets differing in surface gravity, atmospheric compo- 
sition, metallicity, cloud coverage and surface type. We do not 
consider polarized light in this paper and leave such a study for 
future work. Our primary goal is not to refine the instrument c on- 
cept given in Boccaletti et al. (2012). Firstly, we investigate the 
exoplanet detection space realistically, and then we may use our 
results to update the design. Our work will also be beneficial to 
other space coronagraph proposals. In Sect. |2] we describe the 
numerical model and the assumptions we use for our study. In 
Sect. |3] we analyze the effects of speckle noise, read-out noise, 
exozodi and photon noise on the performance. We specify what 
kinds of planets can be detected in each case. We then study the 
impact of the spectral type of the host star In Sect.|4] we define a 
criterion that gives the required flux accuracy that the instrument 
has to produce in order to disentangle spectra for similar planet 
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types (Jupiters, Neptunes or super-Earths), but with different val- 
ues of physical parameters (composition, metallicity or cloud 
and surface coverage). We conclude that our instrument fulfills 
these constraints for most of the considered cases. In Sect.|5] we 
generalize our results to other stellar types to define a parameter 
space of potential targets. 

2. Models 

2.1. Numerical model of the instrument 

The instrument concept of SPICES is designed to provide po- 
larimetric and sp ectroscopic measureme nts in the 0.45-0.90 yum 
range (Fig. fTland lBoccaletti et al.l (l2012h '). To achieve high con- 
trasts required to exoplanet characterization, SPICES combines 
a high-quality off-axis telescope, high-accuracy wavefront con- 
trol, a broad-band coronagraph with small inner working angle 
(IWA; the angular separation of 50% throughput), and optical 
elements and detectors to collect the polarimetric and spectro- 
scopic information. As for the telescope, very low wavefront 
aberrations have already been demonstrated for GAIA primary 
mirrors (~ 8 nm rms on surfaceo, and SPICES can benefit from 
the same technology. The wavefront control is achieved with the 
combination of a focal plane wavefront sensor and a deformable 
mirror (DM). The wavefront s ensing is ach ieved with a self- 
coherent camera (SCC, Galiche r et al.l 120081) . which is a very 
simple modification in the coronagraph design. The SCC spa- 
tially modulates speckles (residual stellar light) that are recorded 
in the science image (no additional channel) and can accurately 
retrieve the wavefront errors (phase and amplitude) that induce 
these speckles. It then drives a DM to coiTect for them. As the 
correction is never perfect, the speckle noise is strongly atten- 
uated, but residual speckles still remain in the image. The SCC 
provides a means to calibrate them and extract the companion or 
disk information, without prior information on the spectrum of 
the latter (Baudoz et al. 2006; Galicher et al. 2010). The DM is a 
64x64 actuator miiTor, as a larger number of actuators incre ases 
the field of view which can be corrected. A vortex coronagraph 
in the focal plane applies an azimuthal phase ram p to the cor- 
rected incident wavefront to cancel the starlight (Maw et et alJ 
l2005h . The voitex coronagraph can be made achromatic over 
a wide spectral bandwidth (-50%, iMawet et al. 2010 ). This 
type o f coronagraph was successfully used by iSerabvn et al.l 
( 1201 Ol) to re-image the HR8799 multi-planet system. Finally, 
an integ ral field s pectrograph (IFS) based on a micro-lens ar- 
ray (An tichi et al.l l2009) allows recording of a dispersed image 
of the corrected field of view. Once a (x, y. A) data cube is recon- 
structed from the detector image, the SCC speckle calibration is 
applied to every spectral channel separately, therefore drastically 
reducing the chromaticity of this device (IGaliche r et al. 2010). In 
the current design, the beam is divided into two branches, each 
assigned to a state of polarization and half of the spectral band. 
An achromatic modulator, at the very beginning of the instru- 
ment, selects the linear polarization direction on the sky that is 
analyzed by the polarizer in each branch. 

We built a numerical code, written in the Interactive Data 
Language (IDL), to model the instrument concept and to simu- 
late the SPICES performance. Our code operates in three steps: 

- Step 1: The simulation of non-coronagraphic and corona- 
graphic image cubes of on-axis (the star) and off-axis (the 
planet(s)) sources. The third dimension of the cubes repre- 
sents the spectral channels. We assume that the image cube 
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Fig. 1. SPICES conceptual baseline. 



extraction from the IFS data is perfect, but correctly account 
for flat field impact at step 2 (Sect. l272l i. 

- Step 2: The normalization of the star and planet spectra, and 
introduction of noise (photon noise, zodi, exozodi, read-out 
noise, flat field). 

- Step 3: The SCC speckle calibration of the images recorded 
at step 2 and the measurement of the planet spectra. 

In Fig. 121 we show an example of simulated images with- 
out detection noise for a single spectral channel (A - 0.675 yum, 
A/l^ 0.013 /urn), after the coronagraph (Fig. |2] left) and after 
the speckle calibration (Fig. |2] middle and right). The image 
field is limited to the zone corrected by the DM (64 x 64 (A/D)^ 
which corresponds to ~6 x 6 arcsec- at A = 0.675 //m). The cor- 
recte d zone size is set by the linear number of DM actua- 
tors (IGalicher et al.ll2010l) . As we introduce both amplitude and 
phase aberrations in the entrance pupil of the instrument and 
SPICES uses a sole DM, the speckles are corrected only in 
one half of the field of view (right half in the images). After 
speckle calibration, the contrast is enhanced as shown in the 
middle and right images. With the current SCC, the calibration is 
more efficient in a diagonal because of the chromatism limitation 
(Galicher et al. 2010). A new version of the technique is under 
study to calibrate speckles in all directions. In the left and middle 
images of Fig.|2] we added two jovian planets of contrasts ~10"^ 
and ~10"^ at 2.2 and 5.4 A/D (red circles). The closest planet 
is detected in the coronagraphic image (left) with a SNR^25. 
The farthest planet can be seen in this image if its position is 
known. However, it is not possible to claim a detection in this 
image (SNR ^ 5), whereas it is detected with a SNR ^ 600 after 
calibration (middle). 

2.2. Instrument assumptions 

In terms of contrast, SPICES has to reach values as low as 10"^ 
at 2 A/D and 10"'° at 4 A/D (Boccaletti et al. 2012) to pro- 
duce interesting science results. Such performance is achieved 
assuming the parameters and the requirements on noise that 
are listed in Table [1] We consider these values in the numer- 
ical simulation and discuss some of them in this section. The 
star is assumed to be perfectly centered onto the coronagraph 
focal plane mask. Since the pointing accuracy is a critical as- 
pect, SPICES will include a dedicated procedure for a precise 
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Fig. 2. Central part of images without detection noise produced by the simulation, after the coronagraph (left) and after the speckle 
calibration by the self-coherent camera (middle and right). The image size is ~64 x 64 (A/D)^ {A - 0.675 fj.m). In the two left images, 
there are two jovian planets of contrasts ~10"^ and ~10"^ (red circles). In the right image, we indicate the calculation area of the 
profiles shown in Sect.|3]with a white dotted line. The intensity scales are linear and identical. 



Table 1. Assumptions used for the instrument simulations. 



Parameter 



Value 



Telescope diameter 

Spectral bandwidth 

Spectral resolution R(A = 0.675 fim) 

Number of spectral channels 

A/i per channel 

AfD sampling at 0.45 fim 

Wavefront errors (A = 0.675 yum) 

Amplitude aberrations 

Number of actuators on the DM 

Phase estimation by the SCC 

Throughput 

Quantum efficiency 

Maximum integration time 

Read-out noise 

Full well capacity 

Time per single exposure 

Flat field accuracy 

Zodiacal light 

Exozodi level 



1.5 m 

0.45-0.90 //m 

50 

35 

-0.013 yum 

3 pixels 

15 nm rms (f-^ PSD" ) 

0.1% (f"' PSD") 

64x64 

perfect 

23% 

70% 

200 h 

0.2 e- rms per pixel 

300 000 e- 

1000 s 

0.5% 

V = 23. 1 mag arcsec"^ 

1 zodi 



Notes. The parameter values are discussed in Sects. 12711 12. 2l and l3. 21 
*"' f is the spatial frequency of the optical defect, PSD refers to power 
spectral density. 



control of the coronag raph alignment at the level of ~0.2 mas 
dBoccaletti et al.ll2012l) . Our simulations showed that this value 
keeps the speckle noise at a level of ~10 '° between 2 and 
4 A/D, which is below the photon noise (Sect. 13.31 1. In the op- 
tical design of the instrument, the whole spectral band is split 
into two branches, the first branch covering the 0.45-0.7 fim 
band (AA/A ^ 43%) and the second branch the 0.65-0.9 jim band 
{AA/A^32 %). The overlapping is for calibration purposes. We 
assume that all phase and amplitude abeiTations are located in 
planes conjugated to the instrument pupil and we use the ma- 
trix direct Fourier transform (Soummer et al. 2007) to propa- 
gate the light. Fresnel p ropagation e ffect will be included us- 
ing the PROPER Ubrary (lKristll2007h . The main impact will be 
a partial modificati on of the speckle pattern with wavelength 
('Marois etaT]l2006l) . This will not impact the SCC wavefront 
estimation since each spectral channel is treated separately, but 
it will red uce the efliciency of the wavefront correction with a 
sole DM dShaklan & Greenll2006l) . A sUghtly undersized Lyot 



stop is used for the coronagraph (95% of the pupil diameter). 
We consider that the vortex coronagraph performance is achro- 
matic for SPICES' spectral bands. CuiTent vortex coronagraphs 
are limited to contrasts of ~4 x 10"** for a 20% bandwidth, but 
strong eff'orts have been made in the past few years to develop 
them in laborator y, test them on the s ky and further improve 
their performance (ISerabvn et al.ll201 ll) . 15 monochromatic im- 
ages are co-added to simulate each spectral channel image. The 
SCC requires an oversampling with respect to Shannon's crite- 
rion. The pixel number per spatial resolution element (A/D) is 3 
at the minimum wavelength. The SCC fringes are tilted at 45° 
with respect to the pixel grid. We use a f"^ power law (f is the 
spatial frequency of the optical defect) for the power spectral 
density (PSD) of the phase aberrations, since it quite well re- 

I Toduces the aberratio ns measured on the VLT and HST miiTors 
Borde & Traubl l200OT. Simulations indicate that the amplitude 
aberrations will have to be <0.1% in order to meet SPICES' re- 
quirements if they follow a f ' PSD, but this value w ill be more 
stringent if the law exponent is >1 ( Galichetjl2009l chap. IV.2, 
Fig. IV.2. 2). SPICES' optical aberrations are expected to evolve 
very slowly with time. SPICES will be located at the L2 point, 
which is believed to be a very stable environment. This may be 
confirmed by the GAIA and JWST missions. Assuming such a 
stable environment, we plan to allocate a significant amount of 
time at the beginning of the mission for the purpose of accurately 
estimating SPICES' aberrations. Then, the DM will compensate 
for the slow variations. In our simulation, the SCC perfectly 
estimates for the wavefront abeiTations (phase and amplitude). 
The perfectly estimated wavefront is projected onto the 64 x 64 
DM using the method of energy minimization in the pupil plane 
(Horde & Traub 2006). Th e DM influen ce functions are modeled 
by adapting the formula of i Huang et al. (2008) to fit the param- 
eters of a realistic DM. Finally, the numerical noise introduced 
by the extraction of individual spectra from the IFS to build data 
cubes is assumed negligible. 

In st ep 2, blackbody spe ctra for the star a nd planetary spec- 
tra froml Cahov et al.l (1201 Ol) and'Staml (l2008b for the planets are 
introduced. The latter are discussed in detail in Sect. 12.41 Photon 
and read-out noise, flat field variations, and zodiacal and exo- 
zodiacal light are accounted for. The instrument throughput is 
set to 23%, considering ~15 optical surfaces of reflectivity 90% 
from the primary miiTor to the detector, and the quantum ef- 
ficiency of the detector is 70 %. Using a n algorithm to correct 
for cosmic ray contamination, iRobbertd (12009) found that sin- 
gle exposures of 1 000 s will keep cosmic-ray induced glitches 
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negligible with respect to a read-out noise level of a few e- rms 
per pixel. We adopt this single exposure for our simulations and 
we do not account for glitches. The detector flat field is mod- 
eled as gaussian noise with a mean of 1 and a rms of 0.5%. The 
IFS spreads the spectrum of an object point on different detec- 
tor pixels. As a consequence, every pixel of the (x, y, A) cube 
is affected by a specific flat field. We use the measured values 
of zodiacal light from Giavali sco et al. (2002). An exo-zodiacal 
disk with a 60° -inclination with respect to a face-on orbit and 
a 45 "-orientation from the horizontal direction is simulated us- 
ing the Zodipic algorithm ( Kuchner 2004). The inclination value 
is the statistical median assuming a uniformly random orienta- 
tion. The orientation corresponds to the SCC fringe. The actual 
scientific strategy of SPICES requires prior knowledge of the 
orientation of the planet orbit. This has to be achieved by com- 
bining astrometric measurements with RV data. The present and 
near-future instruments like VLT/PRIMA and GAIA (we note 
that stars with V < 6 are too bright for the latter) can in principle 
provide the information. We also plan to use a new version of 
the SCC that could enlarge the high-contrast part of the image 
and thus, could relax the constraints on the orbital knowledge. 

2.3. Selecting spectral bandwidth and resolution 

In this section we discuss the different parameters we have to 
take into account for the choice of SPICES' spectral bandwidth 
and resolution. To determine the spectral bandwidth, we first 
examine theoretical spectra representative of Jupiter, Neptune 
and terrestrial atm ospheres (Fig. |3] giant planet models from 
ICahov et all (12010' ) and Earth models from IStamI (12008.) ). The 
spectra are calcula ted for visible w avelengths: .35-1 fim for 
ICahov et all ilQKJti and 0.3-1 ^m for lStaml(i2008l) . We note that 
a 0.45-0.90 fim bandwidth ofifers a good compromise and en- 
ables to measure Rayleigh scattering at the blue wavelengths as 
well as the molecular absorption bands in the red part. The trade- 
off" for the long-wavelength cut-off" results from a technological 
limitation, since visible detectors have weak efficiencies above 
0.90 fim. Therefore, SPICES will not measure the wide and 
strong water absorption band at 0.94 fim in the spectrum of Earth 
analogs (Fig. [3]l. Although deep absorption bands are easily 
identifiable in model spectra, their depth would be hard to mea- 
sure given the SPICES performance. However, other but shal- 
lower water bands at 0.72 and 0.82 /im are also present in plan- 
etary spectra. Similarly, giant planets feature a strong methane 
band at 0.89 //m at the boundary of SPICES' bandpass, but 
there are weaker bands at, for instance, 0.62, 0.73 and 0.79 /im. 
Measuring absorption bands at different wavelengths allows to 
infer the gas abundances, if the cloud top altitudes can be de- 
rived froina known gas which is well mixed in the atmosphere 
(IStamI l2008l and references therein). Polarimetry combined 
with flux could also help to break the degeneracy ('Sta m et aTl 
|2004 . For the telluric planets, Des Marais et al. (2002) define 
the spectral bandwidths and list the molecules that exoplanet 
missions should address: molecular oxygen (O2), ozone (O3), 
water (H2O), methane (CH4) and carbon dioxide (CO2). The 
spectral range of SPICES (0.45-0.90 fim) permits to measure all 
these molec ules ex cept CO2. In addition to atmospheric gases, 
ISeager et al.l(l2005h have emphasized the scientific interest of the 
detection of surface features like the "red edge" (the rise of the 
clear Earth spectrum beyond 0.7 yum in Fig. |3]l. Nevertheless, 
several studies showed that this measurement is difficult for 
the Earth itself (e.g., Woolfetal. 2002; Arnold etal. 2002; 
iMontafies-Rodriguez etaill2005l) . ISeager et al.1 (l2005l) note that 
the "red edge" should be detected with molecular oxygen to be 
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Fig. 3. H igh-resolution albedo spectra of Jupiter, Neptune 
(ICahov e t al. 2010) and two Earth analogs (Stam 2008). The 
curves of the Jupiter and Neptune are vertically shifted with re- 
spect to the actual albedo for the sake of clarity. The dashed hor- 
izontal lines indicate the reference position for the Jupiter spec- 
trum (black) and the Neptune spectrum (red). 



sure that it is related to vegetation, because minerals may present 
a similar feature but at different wavelengths. From the compar- 
ison of widths of SPICES' spectral channels and of molecular 
bands, we find that a spectral resolution of at least 50 is required 
to identify the main bands of the spectra of giant planets as well 
as super-Earths (Tables |3] and El and Figs. |4] and 5). Our analysis 
confirms previous results ([Schneider et al.ll2009b . 



2.4. Planetary atmosphere models 



The following sections 12.4.11 and 12.4.21 describe the models we 
use in our simulation. In this paper, we consider two sets of mod- 
els, one for giant planets and one for telluric planets. Cahoy et alj 
(2010)'s models cover a large range of star-planet separations 
and planet metallicities. IStamI (l2008l) 's models assume different 
planet surfaces and atmospheric structures (cloudy and clear at- 
mospheres). Table |2] summarizes the main parameters we con- 
sider for these models. We set a maximum mass of 10 Earth 
masses (Me) for the super-Earths, which corresponds to a maxi- 
mu m radius of 2 . 5 Eart h radii (Re) from the mass-radius relation 
of iGrasset et al.l (l2009h . Recall that Stam's atmosphere models 
assume an Earth-like planet. 

From the degraded spectra of Figs. E]and|5](see Sect. l2.3l for 
the discussion of the spectral resolution), the planet/star contrast 
is expected to be ~ 10^^-10"'° for Jupiters and ~ 10^^-10 " for 
Neptunes and super-Earths, depending on the separation from 
the star and on the wavelength. Note that Cahov et al. (201(|) 
and lStamI (l2008h present geometric albedo spectra, while we plot 
contrast curves using the following formula: 



^2 



(1) 



where A is the wavelength, C(A) is the planet/star contrast spec- 
trum, A{A, a) is the planet albedo spectrum at phase angle a, 
Rp is the planet radius and a is the star-planet separation. Phase 
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Table 2. Parameters of the exoplanetary atmosphere models used in this paper. 



Planet 



Separations (AU) Radius Metallicity (solar units) Atmospheric structure 



Surface type 



Jupiter 


0.8, 2, 5 and 10 


iRj 


lands 


Neptune 


0.8, 2, 5 and 10 


IRn 


10 and 30 


Super-Earth 


1 


2.5 Re 


- 



0, 50 and 100% clouds Forest, ocean and forest-ocean mix 



Notes. Rj, Rn and Rg refer to Jupiter, Neptune and Earth radii respectively. 
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Fig. 4. Theoretical models of Jupiter-Uke (left p anel) and Neptune-like (right panel) atmospheres degraded to the resolution of 
SPICES (models adapted from lCahov et al.ll2010h . 



angle is the angle star-planet-observer: when the planet is be- 
hind the star a = 0°, and when the planet is in front of the star 
a - 180°. In this work, we only consider the case of a planet at 
its maximum elongation from the star, a - 90°. The flux reflected 
by a planet depends on both its albedo and its radius (Eq. ([TJ). 
To determine the albedo from observations, we need to inde- 
pendently estimate for the radius. We can u se theoretical mass- 
radius relations (e g., iFortnev et all l2007l iBaraffe et all l2008t 
iGrasset et al.l l2009h . The polarimetric capabilities of SPICES 
could also help to establish the planetary properties without prior 
information about the planet's radius (Stametal. 2004; Sta^ 
[2008). This will be subject for future work. When we will study 
SPICES' ability to retrieve the planet properties from measured 
spectra in Sect. H] we assume that the radius is known. 



Table 3. List of molecules observable at R ■ 
and Neptune models of Fig.|4] 



: 50 in the Jupiter 



Molecule 


Appr 


aximate A 


(jj.m) 


Note 


CH4 




0.54 






CH4 




0.62 






H2O 




0.65 




for 0.8-AU models 


CH4 




0.66 






CH4 




0.73 






CH4 




0.79 






H2O 




0.82 




for 0.8-AU models 


CH4 




0.84 






CH4 




0.86 







2.4.1 . Theoretical spectra of Jupiter and Neptune analogs 



We use the models of ICahov et al.l (l20IOl) . who calculate atmo- 
spheiic structures of old (~4.5 Gyr) Jupiter and Neptune analogs 
in radiative equilibrium with the radiation of a solar-type host 
star at separations of 0.8, 2, 5 and 10 AU, for different metal- 
licities and for optical wavelengths (0.35 to 1 /im). This range 
matches the separations of the planets that small space coro- 
nagraphs can potentially observe. Table |3] lists the main spec- 
tral bands observable with a spectral resolution of 50. Several 
theoretical spectra of Jupiter and Neptune analogs are shown 
in Fig. H] Unlike young planets, which are dominated by ther- 



mal radiation, the star-planet separation drastically alters the 
structure and composition of mature planetary atmospheres. 
Therefore, a simple scaling of the amount of reflected light with 
distance is not sufficient to model realistic spectra and to derive 
the actual performance of a mission like SPICES. If planets are 
too warm for any molecules to condense into clouds, their spec- 
tra are dominated by Rayleigh scattering. This is illustrated for 
the case of a separation of 0.8 AU in Fig. |4] At 2 AU, blight wa- 
ter clouds form and dominate the atmospheric opacity all over 
the spectrum. At 5 AU, ammonia clouds form above the water 
clouds. At 10 AU, the same clouds form but at a deeper pres- 
sure level, and Rayleigh scattering again dominates the reflected 
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Fi g. 5. Same a s Fig. |4]but for the terrestrial atmosphere models 
of [staml(l2008h . 

Table 4. Same as Table[3]but for the terrestrial planet models of 
Fig.S 



Molecule 


Approximate A 


(jum) 


Note 


O3 


0.5 - 0.7 




the Chappuis band 


O2 


0.69 




B-band 


H2O 


0.72 






O2 


0.76 




A-band 


H2O 


0.82 







Table 5. Star-planet separations for giant planets extrapolated to 
other stellar types assuming Eq. ^. The values for the solar-type 
star are those modeled by Cahoy et al. (2010). 



Spectral type 


Luminosity (L©) 


Separations (AU) 


AO 
FO 


28 
4.8 


4.2 10.6 
1.8 4.4 


26.5 53.0 
11.0 21.9 


■ G2 1 


^^^■^^^P 


0.8 2 


5 10 


KO 
MO 


0.45 
0.09 


0.5 1.3 
0.24 0.6 


3.4 6.7 

1.5 3 



water cloud layer, and surfaces completely covered by either for- 
est or black ocean with a Fresnel reflecting interface (Fig. |5]l. 
The water clouds are optically thick and located in the tropo- 
sphere. The atmospheric absorbers are water, molecular oxygen 
and ozone (Table HJ. The model albedo of vegetation presents 
two main features: a local maximum between 0.5 and 0.6 jum, 
which is due to two absorption bands of chlorophyll at 0.45 and 
0.67 yum, and the increase of the albedo beyond 0.7 fj.m due to the 
"red edge" (Sect. I2.3l l. The appearance of the "red edge" in the 
planet spectra depends strongly on the cloud thickness and cov- 
erage, but it still alters the spectrum shape for a partial cloud cov- 
erage despite their large optical thickness (Fig. |5] see the spec- 
tra of the 50% cloudy ocean and 50% cloudy forest planets in 
blue and green dashed lines respectively). At short wavelengths 
Rayleigh scattering dominates while at long wavelengths scatter- 
ing by clouds is the most important process. 'Stam| (l2008b uses a 
weighted sum of homogeneous models to simulate a quasi hori- 
zontally inhomogeneous model representative of the Earth (70% 
of the surface covered by ocean and 30% by forest) with dif- 
ferent cloud coverages. In this paper, we use different weighted 
sums for simulating three cloud coverages (0, 50 and 100%) and 
three surface types (forest, forest-ocean equally mix and ocean), 
as indicated in Table |2l 



flux at short wavelengths. The planet is ~4 times fainter than the 
Jupiter at 5 AU, as expected in the case where the flux decrease 
follows an inverse square power law of separation. For all cloudy 
planets (2, 5 and 10 AU), the clouds are optically thick. The 
Neptune spectra exhibit the same but stronger absorption bands 
as the Jupiter spectra. Increasing the planet metallicity usually 
decreases its albedo. The exception is the 2-AU Jupiter model 
for short wavelengths because for this case the water c louds are 
high in the atmosphere and thick. ICahov et aP (1201 Oh note that 
the metallicity increase produces larger differences between the 
Jupiter spectra than between the Neptune spectra for separations 
of 5 and 10 AU. Methane bands dominate the spectra over all the 
bandwidth (e.g., 0.62, 0.73, 0.79 and 0.89 //m). Their depths de- 
pend on the nature of the light-scatteiing particles (gases, clouds, 
aerosols). 

2.4.2. Theoretical spectra of telluric planets 

In this section, we summari ze the main features of the telluric 
planet models. IStamI (l2008h uses a pressure-temperature verti- 
cal profile of the Earth to derive flux and polarization spectra 
for several surface and cloud coverages and for a spectral range 
between 0.3 and 1 yum. We recall that we do not consider polar- 
ization effects in this paper and that these mod els are releva nt for 
a separation of 1 AU around a solar- type star. IStanil (|2008|) con- 
siders atmospheres with and without a horizontal homogeneous 



2.4.3. Planet contrasts as a function of the stellar type 

The atmospheric structures of the models were calculated by as- 
suming ~4.5-Gyrplanets in radiative equilibrium with the flux of 
a Sun-like parent star. These models can be transposed to other 
stellar types assuming flux conservation with the standard for- 
mula: 



AnRlaT% = 



(l-AB);r/?2 * 
' 47rfl^ 



(2) 



where Rp is the planet radius, cr is Stefan's constant, Ti,q the 
planet atmosphere equilibrium temperature, Ag the planet Bond 
albedo and L* the host star luminosity. This formula does not 
account for the effects of the wavelength dependerice of the star 
emission on the atmosphere (iMarlev et al.lll999t iFortnev et al.l 
2007). 

Because we are using models with discrete values (especially 
separations and stellar luminosity), we cannot extrapolate the 
planet spectra to any separations around any stars. Instead, we 
calculate the correspondence between separations and stellar lu- 
minosity, considering T^q and Ab only depend on the incident 
stellar flux at the planet position L*/a^. Therefore, Eq. (13) be- 
comes: 



Li, cc a 



as,. 



ac2 



(3) 



where as is the star-planet separation for a star of spectral type 
S p, flc2 the star-planet separation for a G2 star and Ls , the star 
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Fig. 6. 5-cr detection profiles of the instrument contrast achieved 
by SPICES for all spectral channels, before (blue solid lines) and 
after (red dashed lines) the SCC speckle calibration. 



luminosity in solar units. For example, a Jupiter at 2 AU from 
a solar-like star would have the same atmospheric structure as 
a Jupiter at ~10.5 AU from an AO star. Table|5]gives the corre- 
spondences for different star-planet separations and stellar types. 
Substituting Eq. Q into Eq. ([TJ, we obtain the contrast of a 
planet around a host star of type S j,: 



C(A) = A(A, a) 



Ri 



a 02 ^s„ 



(4) 



We use Eqs. (O and (|4|i to derive the star-planet separations and 
contrasts in Sect. [3] While this calculation is a fairly good es- 
timation for the giant and cloudy telluric planet models, it is 
less accurate for the clear terrestrial mod els. For stars cooler 
than the Sun. lWolstencroft & Raved (|2002|) suggest that the "red 
edge" could indeed be shifted towards wavelengths redder than 
1 //m, if the photon number involved in the photosynthesis pro- 
cesses is greater than for the mechanism operating on Earth. 
If this hypothesis is verified, the "red edge" will be outside 
the bandwidth covered by SPICES and undetectable. However, 
[Kiang et al. (2007) warn that a theory predicting the "red edge" 
wavelength for a given stellar type assuming the same mecha- 
nism as on Earth is still missing. 



3. Performance in detection 

The instrument model presented in Sect. 12.1 [ provides an estima- 
tion of the achievable contrast map in the field of view for each 
spectral channel. In this section, we estimate the average con- 
trast that is reached in the darkest area of the field of view (area 
contained by the dotted lines in Fig.|2]l as a function of the angu- 
lar separation from the central star. As we explained in Sect. 12.21 
we assume that the planet orbital parameters are known and its 
position can be matched with the orientation of the corrected 
area. 



Fig. 7. 5-cr detection profiles of the instrument contrast averaged 
over the SPICES bandpass (solid lines) and compared to aver- 
aged planet contrasts calculated for a solar-type star (symbols). 
For comparison, we also plot the mean contrast of the zodiacal 
light (horizontal dashed line) and a 1-zodi exodisk (dot-dashed 
curve), both for a star at 10 pc. 



3.1. Impact of speckle noise 

Figure |6] shows SPICES' contrast curves against angular sepa- 
ration (A() - 0.675 fim) for all spectral channels, before and af- 
ter the speckle calibration by the SCC. We see that the speckle 
subtraction is necessary to reach the requirement of ~10"'° at 
a few A/D (Sect. lZSl l. The wavelength dispersion of the perfor- 
mance is due to the phase abeiTation dependence on wavelength 
(oc/l~^) and the SCC calibration dependence on spectral reso- 
lution (iGaUcher et al.l l2010). Recall that we set the same band- 
width for all channels so spectral resolution increases with wave- 
length. The steep increase of the detection limit around 32 A/D 
coiTesponds to the DM cut-off spatial frequency. This cut-off and 
the more efficient SCC calibration at small separations explain 
the degradation with angular separation. 

To have a clear and simple view of the performance, we plot 
radial profiles averaged on all the spectral channels against the 
physical separation in AU for three star distances in Fig. |7] The 
performance is limited at short separations by the coronagraph 
IWA {~2 A/D, Sect. 12.1b . and at large separations by the size of 
the DM corrected area (cut-off at 32 A/D). As we express the 
X-axis in AU, the contrast curve scales with the star distance. 
We also plot the positions of Jupiter of 1 solar metallicity and 

1 Jupiter radius (Rj), Neptune of 10 solar metallicities and 1 
Neptune radius (Rn), and 2.5-Re cloudy planets. We assume that 
the super-Earth properties do n ot evo l ve wit h separation for the 
considered range. We recall that lStanJfcOOSh uses a temperature- 
pressure profile of the Earth so the model is relevant for a sep- 
aration of 1 AU. For stars at 20 pc, the farthest Jupiters (5 and 
10 AU) and Neptunes (5 AU) are detected with SNR>5. The 
10- AU Neptune is below the curve (SNR ^ 4.5). For a 10-pc star, 
the same planets are still detected as well as planets as close as 

2 AU. For the closest star we consider (3 pc), planets as close as 
1 AU are very well detected and characterized (Sect. 2). In this 
case, the 10-AU Jupiter and 5-AU Neptune are not detected with 
the current instrument design because of the speckle calibration 
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Fig. 8. 5-cr detection profiles accounting for the photon noise of Fig. 9. 5-cr detection profiles for a Sun analog at 3 and 10 pc and 
G2 star at 3 (blue lines) and 10 pc (red lines), and for read-out exo-zodiacal intensities ranging from 1 to 1 000 zodis, assuming 
noise levels from to 3 e- rms per pixel. a 0.2 e-/pixel rms read-out noise and photon noise. 



degradation with angular separation. Solutions may exist to im- 
prove the detection at these larg e angular separations (>18 AjD) 
such as spectral deconvolution (ISparks & Fordll20()2l) . 

In addition, we plot the averaged contrast levels for zodia- 
cal light and a 1-zodi exodisk (star distance of 10 pc). Recall 
that they increase as the square of the star distance. Although 
the instrument concept can reduce speckle noise below a con- 
trast of 10~'°, the final performance is limited by these extended 
background sources. The considered exo-zodiacal disk limita- 
tion is at the level of Neptunes and super-Earths. Therefore, to 
obtain a correct estimation of planet fluxes, the zodiacal and exo- 
zodiacal contributions must be carefully calibrated and removed. 
The exo-zodiacal disk intensity has been identified as critical 
for Earth-twin detection. To address this question, exo-disk sur- 
veys have been recommended to prepare target lists of faint 
exo-disks for space nulling interferometers like the Terrestrial 
Plan et Finder Interferometer an d Darwin (p^unine et al.l I2008t 
ICou de du Foresto et al. 2010; Ha tzes et al.ll2010l) . In the case of 
SPICES, the problematic is different as we do single-aperture 
imaging. The exact procedure to account for zodiacal and exo- 
zodiacal contributions remains to be defined although it is a data 
reduction issue, which is beyond the scope of this paper In the 
following, we consider the model distribution of both zodiacal 
and exo-zodiacal intensities can be subtracted from the data. The 
photon noise of these contributions may still limit the contrast 
performance (Sect.|32). 



3.2. Impact of detection noise 

Due to a small primary mirror (typically 1.5 m in diameter) to- 
gether with a spectral resolution R - 50, space coronagraphs like 
SPICES will be limited by photon noise (from the stellar back- 
ground and the planet) or read-out noise for most of the tar- 
gets, because very long exposures would be required to reach 
the instrument limitation at ~ 10"'°. For instance, ~10000 h are 
needed to achieve ~10"'" at 5 cr for a G2 star at 10 pc. From 
this exposure time and assuming that the noise follows the pho- 
ton noise behavior (proportional to the square root of the photon 



number), we estimate the total integration time required to de- 
tect a super-Earth at SNR = 5, for typical values of distance and 
contrast accessible to SPICES. It is 200 h for a planet of contrast 
2.5 X 10"'° at 2 AU around a solar-type star at 5 pc. No detailed 
study of the maximum integration time per target has been car- 
ried out to date for SPICES. In this paper, we set the maximum 
integration time per target to 200 h. This is a good trade-off be- 
tween achieving high contrasts and observing a large number of 
targets during the mission. 

We present the impact of photon and read-out noise on 
SPICES' performance for the case of a G2 star at 3 and 10 pc in 
Fig. |8] We note that the read-out noise is a major limitation for 
the furthest star but not for the closest. This is due to the fact that 
the full well capacity of the detector (Table[T]i is not filled after a 
1 000-s exposure in both cases. The number of single exposures 
and read-out noise level are thus the same, but the photon count 
is greater for the closest star We note that the dozen of stars lo- 
cated within 3 pc have types later than G (the exception being 
Sirius). Thus, the read-out noise will not be a fundamental limi- 
tation for close stars. We base the read-out noise requirement on 
the farthest star and set its value to 0.2 e- rms per pixel. Electron 
multiplying CCDs can achieve such a low read-out noise, and 
a large set of devi ces have been qua lified for space during the 
GAIA preparation (ISmith et al.ll2006h . 

As indicated in the previous section, an exo-disk can pre- 
vent the detection of faint planets if its photon noise becomes 
too important. To help to prepare a target list, we estimate the 
exozodi level that may hamper the detection of SPICES' targets. 
Figure |9] presents the performance for different exozodi levels 
and two distances of a solar-type star. The read-out noise is set 
to 0.2 e- rms per pixel. We assume that the exo-disks have no 
structure and can be subtracted out from the data to the precision 
imposed by photon noise. We find that the exo-disk photon noise 
does not significantly limit the performance up to 10 zodis, but 
begins preventing the Neptune and super-Earth detection when 
larger than 100 zodis. For exposures shorter than 200 h (more 
dominant photon noise), the acceptable exozodi level is lower. 
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Fig. 10. 5-0" detection profiles for a 200-h exposure compared 
to averaged planet contrasts for G2 (top), AO (middle) and MO 
(bottom) type stars and several star distances. 



3.3. Impact of stellar type 

Assuming a generic exposure time of 200 hours and an exo- 
zodiacal disk of 1 zodi, we test the detectability of planets for 
several stellar types. Figure [10] presents the 5-cr detection lev- 



els for three stars of type G2, AO and MO and several dis- 
tances (3, 10, and 20 pc). We consider the star-planet separa- 
tions given in Table|5](they scaled by L^ ), and we apply Eq. (HI 
to find the wavelength-averaged planet contrasts. For a G2 star 
(Fig.[TOl top), SPICES can potentially access jovian planets, icy 
giants and super-Earths with separations smaller than ~6, ~3 and 
~2 AU respectively. No planet is detected at distances larger than 
~10 pc. The detectability of the closest planets (0.8 to 2 AU) 
will be limited by the coronagraph IWA (~2 A/D) for stars at 

4 to 10 pc. As for an AO parent star (Fig. [10] middle). Jupiter- 
like planets are not detected when further than 10 AU for stars 
within 20 pc. Th e separation ranges considered by Cahov et aP 
(1201 Oh and lStamI feOOS.) in the case of a G2 star do not allow to 
study planets at separations shorter than ~4 AU around an AO 
star. Nonetheless, we can roughly estimate that Neptune analogs 
and super-Earths can be detected in the range 2-4 AU for stars 
closer than 10 pc. There are no A stars within 5 pc (except for the 
Sirius binary system) so we do not plot the 3-pc curve. Finally, 
the very close M stars (Fig. [10] bottom) at 3-5 pc are of great 
interest for detecting Jupiter-like planets in the 0.5-4 AU range 
as well as super-Eaiths in close orbits at 0.5-1 AU. 

From the current exoplanet database (ISchneider et al.lUoi ll) . 
we assess that only a handful of known extrasolar planets match 
the limitations described here. However, we note that RV surveys 
are not complet e in the c ase of early and late stellar types (A and 
M in particular. lUdrv & Santos 2007i) and at long periods (a few 
AUs) even for nearby stars. 

4. Performance in spectrometry 

In this section, we analyze the SPICES performance more in de- 
tail, by estimating SNRs of the measured planetary spectra. The 
objective is to set the constraints on the SNR to allow the dif- 
ferentiation between planetary models: impact of physical star- 
planet separation and metallicity for the Jupiter and Neptune 
analogs (Sects. I4~2l and l4.3l l, and variations due to cloud and sur- 
face coverage for the rocky planets (Sect. l4~4t . 

4.1. Criterion of characterization 

We first explain our characterization criterion. The underlying 
question is to know if a measured spectrum S is reproduced by 
one of two model spectra noted M, (/ refering to the model index, 
;= 1,2). This depends on the noise of the measured spectrum 
A^ - 5 /SNR. We define the following criteria of comparison of 

5 to Mr- 



criti - median^ 



S(A) - MjjA) 
N(A) 



(5) 



where the median is calculated over the spectral channels. We 
choose the median because it accounts well for the overall qual- 
ity of a spectrum. For a given measurement (S , A^), the model 
which best matches the measured spectrum gives the lowest cri- 
terion value. Substituting the definition of A^ to express cr/f, as 
a function of SNR and assuming the latter is nearly constant on 
the measured spectrum, we can write: 



criti = median^ | | x SNR 

S{A) 



(6) 



In this paper, we use theoretical models to assess the instrument 
performance and to set the SNR required to measure differences 
between them (S = Mi and M, = M2). We consider that the mod- 
els are differentiated when their difference is 10 times above the 
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Table 6. Values of SNR,- derived from Eq. ©. 



Planet 


Parameter 


SNR, 


Note 


Jupiter 


0.8/2 AU 


15 




Jupiter 0.8 AU 


metallicity l/3x 


30 




Jupiter 2 AU 


metallicity l/3x 


30 


CH4 bands 


Jupiter 5 AU 


metallicity l/3x 


30 


CH4 bands 


Neptune 


0.8/2 AU 


15 




Neptune 0.8 AU 


metallicity 10/30x 


30 




Neptune 2 AU 


metallicity 10/30x 


25 


CH4 bands 


Forest Earth 


0/50/100% clouds 


25 


blue channels 


Ocean Earth 


0/50/100% clouds 


25 


blue channels 


Clear Earth 


0/50/100% forests 


12 


red channels 


50% cloudy Earth 


0/50/100% forests 


30 


red channels 


Cloudy Earth 


0/50/100% forests 


220 


red channels 



Table 7. Maximum star distance at which SPICES resolves the 
planet separation at the central wavelength of the bandwidth. 



Planet separation (AU) Star distance (pc) 



0.8 
1 

2 
5 
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Fig. 11. Evolution of the median value of SNRm (see text) as 
a function of exposure time for Jupiter planets (symbols). The 
curves are power-law fits. We plot 1-cr eiTorbars. 



measured noise so cr/f, = 10. Inverting Eq. (|6]l and calling SNRr 
the value of the required SNR, we finally obtain: 



SNRr = 10 X 



1 



median^ ( 






(7) 



As an example, when we study the metallicity (Sects. 14.21 and 
I4.3l l, M\ will be the lowest metallicity model and M^ the high- 
est metallicity model. We consider different spectral bandwidths 
to determine SNR, (Eq. (|7])) when analyzing different plane- 
tary properties. For instance, metallicity strenghtens the bands 
(Fig. m, and cloud and surface coverage alter the spectrum at 
blue and red wavelengths respectively (Fig.|5). We give the cal- 
culated SNRi in Tableland we specify the spectral channels in 
the last column. We stress the point that the SNRr values cor- 
respond to the spectrum of the brightest planet considered for 
each analyzed property: for instance, the low-metallicity planets 
for the Jupiters and Neptunes (Sects. I4r2] and l4.3t and the cloudy 
planets for the super-Eaiths (Sect. 14.4) . 

In the remainder of this section, we study the ability of 
SPICES to disentangle planetary models. For each planet sep- 
aration, we consider the distance at which the star-planet system 
is resolved at quadrature (Table |7) and we derive the exposure 
time to achieve the SNR^ values quoted in Table |6] We restrain 
the study to the case of a solar-type star and we assume a max- 
imum exposure time of 200 hours (~8 days). We perform our 
simulations for five independent realizations of speckle pattern. 
We then average our results to minimize the impact of an opti- 
mistic or pessimistic speckle pattern. To save computing time, 
we use the same five speckle patterns for all planet cases, al- 
though we randomly change the photon and read-out noise. We 
assume the planet position to be perfectly known. We integrate 
the planet flux within apertures of diameter I A/D for each spec- 
tral channel. This corresponds to the full width at half maximum 
of the point spread function. The 1-cr eiTor bars shown in the 
plots account for the variation of both speckle and noise realiza- 
tions. 



4.2. Jupiter models 

We consider the models of Jupiter analogs described in 
Sect. 12.4.11 for several separations and metallicities. Figure [TT] 
shows the evolution with exposure time of the median SNR mea- 
sured from the simulated data (SNRn,). We represent each data 
point with its corresponding 1-cr error bar. The SNRm depen- 
dence on exposure time may change from one observed planet 
to another as a function of the planet intensity and location in 
the diffraction pattern of the host star. We fit power-law curves 
since we expect SNRm to be proportional to the square root of 
the integration time, if the dominant noise is the photon noise, 
or SNRm to be constant, if it is the speckle noise. We find that 
all exponents are close to 0.5, which corresponds to the case 
of photon noise limitation. SNRm rapidly increases with time 
for the brightest Jupiter models (separations of 0.8 and 2 AU), 
while the growth is slower for the faintest model (5 AU). We use 
Fig.[TT]to derive all the exposure times given in this section. 

As long as astrometric measurements are not available, the 
orbital inclination is unknown for non-transiting planets detected 
by RV. We recall that we use an inclination of 60° which is the 
statistical median for random orientations (Sect. 12. 2K Although 
direct imaging is able to put constraints on this parameter, it 
requires several images in practice covering the orbit and high 
SNRs. The fact that planets are easier to observe at quadrature is 
also not particularly favorable for a precise determination of in- 
clination and then mass. Therefore, in the case of a single obser- 
vation, the physical separation of the planet to the star is poorly 
const rained, particularly sin ce eccentricity could also be signif- 
icant (lUdrv & Santosll2007h . We can thus confuse a giant planet 
close to its star with a large planet at large separation if their 
projected separation is the same on the image. In addition, for 
an eccentric orbit, a planet's albedo can depend strongly on the 
orbital position: the planet can be almost cloud-free near peri- 
helion and covered by clouds near aphelion. The spectroscopic 
characterization could help to break degeneracies in these pa- 
rameters, if the spectral differences are large enough to be de- 
tected. Considering the theoretical models, giant planet spectra 
mostly diff'er in the blue, where Rayleigh scattering dominates 
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Fig. 12. Spectra of Jupiters at 0.8 and 2 AU as they would be 
measured by SPICES (symbols) and the corresponding model 
spectra (curves). 



for a planet at 0.8 AU. The application of the criterion defined 
in Eq. Q to giant planet spectra indicates that SNRi ^15 per- 
mits to distinguish between the atmospheres of two giants at 0.8 
and 2 AU respectively (Table |6]l. This performance is achieved 
in ~30 min for a distance of 4 pc (Fig. fTTT l. which corresponds 
to the upper limit at which a separation of 0.8 AU is accessi- 
ble to SPICES (Table|2ll. We plot the two spectra as they would 
be measured by the instrument with 1-cr error bars as well as 
the corresponding models in Fig. [12] As expected, the blue half 
of the bandwidth is the region where the two spectra can be 
distinguished with no ambiguity. SPICES will be able to mea- 
sure Rayleigh scattering and estimate the star-planet separation. 
However, we note that these measurements would be possible 
for a few stars only because of the small angular resolution. 

We now study SPICES' ability to measure the spectral differ- 
ences between Jupiter models with 1 and 3 times the solar metal- 
licity for star-planet separations of 0.8, 2 and 5 AU at levels an 
order of magnitude above the noise. As indicated in Table|6] this 
requires SNRr ^ 30 for all separations. For the 0.8-AU Jupiters, 
this value is achieved within the distance for which our instru- 
ment can resolve such a planet (<4 pc). In particular, an inte- 
gration time of 2 h satisfies the criterion for a star at 4 pc, and 
the differences between the measured spectra are 10 times larger 
than the noise over the blue half of the spectra (Fig. [13] top 
panel). Recall that the error bars shown in the figures are at 1 cr. 
Similarly, the SNR criterion is satisfied for the 2-AU planets at 
the maximum distance of 10 pc for an exposure time of ~30 h. 
Metallicity effects are larger in the methane bands for this case 
and we focus on the 0.73-jum methane band and the blue edge 
of the 0.89-//m deep band to disentangle the spectra (Fig. [13] 
middle panel). For a 5-AU Jupiter, the two metallicity cases are 
distinguished in 200 h at a distance of 7 pc in the methane bands 
at 0.62 and 0.73 yum and on the blue edge of the deep feature 
at 0.89 fim (Fig. [13] bottom panel). We note that for the 0.8- 
and 5-AU cases, the metallicity effects can mimic a radius vari- 
ation by shifting the whole spectrum, while for the 2-AU case 
they alter the flux specifically in the absorption bands. We con- 
clude that for resolved systems, SPICES will be able to analyze 
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Fig. 13. Measured (symbols) and theoretical (lines) spectra of 
the 0.8-AU (top), 2-AU (middle) and 5-AU (bottom) Jupiter 
models for 1 and 3 times the solar metallicity. The planet/star 
contrast scale is identical to Fig. [T2] for comparison except for 
the bottom panel. 



metallicity enhancements as small as a factor of 3 for all 0.8- and 
2-AU Jupiter targets around solar-type stars. As planets at 5 AU 
are fainter, they will be accessible only for G2 stars within 7 pc, 
considering the maximum exposure time of ~200 h. 
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Fig. 14. Same as Fig.[TT]but for the Neptune models of Sect. 14. 3 
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Fig. 15. Same as Fig.[T2]but for the 0.8- and 2-AU Neptunes. 
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Fig. 16. Same as Fig. [13] but for the measurement of metalllc- 
ity effects between 10 and 30 times the solar metalllcity of the 
0.8-AU (top) and 2-AU (bottom) Neptunes. The vertical scale is 
identical to Fig.[T5]for comparison. 



4.3. Neptune models 

Figure fT4lis similar to Fig.[TT]but for the Neptune analogs and is 
used to derive the integration times mentioned below. According 
to Table|6l the distinction between the 0.8- and 2-AU planets re- 
quires SNRr ~ 15 in the blue part of the spectral range. Similarly 
to the Jupiter spectra, the spectrum of very close-in Neptunes 
will feature a negative spectral slope, due to Rayleigh scattering, 
which noticeably differs from the nearly flat spectrum of a far- 
ther planet. For the maximum distance of 4 pc at which a 0.8-AU 
Neptune is angularly resolved by SPICES, this value is reached 
in ~50 h. Therefore, the same analysis can be performed for 
closer solar-type targets (as we explained in the previous section, 
Rayleigh scattering would be measurable for a few objects only). 
The simulated measurements for the 4-pc Neptunes at 0.8 and 
2 AU are plotted in Fig. [15] The 0.8- and 2-AU Neptune spectral 
differences can be measured for A < 0.58 //m, whereas the spec- 



tral measurements are degenerate for redder wavelengths given 
the noise level. 

We then test if SPICES can distinguish metalllcity effects 
between Neptunes of 10 and 30 times the solar metalllcity 
(Table HI). W e recall that these values are those studied by 
Cahovetai] (|2010|) . The required SNR, is -30 over the full spec- 
tral range for disentangling the 0.8-AU spectra. It is ~25 in the 
methane bands for the 2-AU planet spectra (Table|6). For the 0.8- 
AU planets, a ~200-h exposure is requested to measure metalllc- 
ity variations for distances as far as 4 pc. The spectral differences 
are detected over the blue channels up to ~0.65 //m (Fig. [16] top 
panel). For a separation of 2 AU, the bottom panel of Fig. [16] 
shows that the metalllcity signatures mainly impact the methane 
bands at 0.62 and 0.66 yum as well as on the edges of the 0.79-/vm 
deep band. We find that SPICES can distinguish the spectra for 
stars within ~6 pc. 
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Fig. 17. Same as Fig. [TT]but for the 2.5-Re planet models dis- 
cussed in Sect. 14.41 
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Fig. 18. Measured (symbols) and model (lines) spectra of 2.5-Re 
forest planets at 1 AU for cloud coverage of 0, 50 and 100%. 
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Fig. 19. Measured spectra of 2.5-Re planets for the clear (top) 
and 50% cloudy (bottom) models with surfaces composed of 
forest, 50% forest-50% ocean and ocean with the correspond- 
ing theoretical spectra. The vertical scale is identical to Fig. [18] 
for comparison. 



4.4. 2.5-Re planet models 

In this section, we analyze SPICES' capability to measure the 
properties of super-Earths, and in particular the effects of cloud 
and surface coverage for planets at 1 AU (their parameters are 
given in Table |2]i. Figure [TT] represents the SNR^ evolution as 
a function of the exposure time for several models of planets 
discussed in this section. The labels "clear" and "cloudy" refer 
to the models with 0% clouds and 100% clouds respectively. The 
star distance is limited to 5 pc (Table |7]i to allow the planet to be 
angularly separated from the star. A few G-type stars would thus 
be accessible to SPICES. 

We first investigate the impact of cloud coverage for both 
forest and ocean surfaces. The influence of clouds is more im- 
portant in the blue (Fig. [18] for surfaces entirely covered with 
forest), because of the strong "red edge" reflection in the red. 



and SNRi ~ 25 (Table |6]l is needed to measure the difference be- 
tween 0, 50 and 100% clouds. At the maximum distance of 5 pc, 
this performance is met in ~150 h for a super-Earth entirely cov- 
ered with forest around a G2 star (Figs. [17] and [18). Systems 
closer than 5 pc will also be accessible for the characterization 
of the cloud coverage. We do not compute the exposure time for 
surfaces covered by ocean, but it will be similar since SNR^ ~ 25 
(Table [6]l and the cloudy spectra are nearly identical for both 
surfaces (green and blue dotted lines in Fig. [5). Figure [TSl also 
shows that the main atmospheric gases (O2 and H2O) are quite 
well retrieved whatever the cloud coverage (SNRn, > 5 and > 15 
respectively). The broad ozone signature is mostly detectable 
when the cloud coverage is large (SNRn, ^ 23). We recall that 
the band wavelengths are given in Table [4] 

Information about the planet surface can be obtained for 
moderate cloud coverage, because it produces noticeable signa- 
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Table 8. Maximum star distance at which SPICES resolves the 
planet separation at the central wavelength of the bandwidth in 
the case of a MO host star. 



Planet separation (AU) 


Star distance (pc) 


0.24 


1.2 


0.3 


1.5 


0.6 


3 


1.5 


7.5 


3 


15 



tures (Fig.[T9ll. The differences between surface types are larger 
for a clear atmosphere, especially in the half red part of the spec- 
tral range. We consider three generic cases of planets with an 
ocean, an equally mixed surface of ocean and forest, and a for- 
est. To separate these cases, SNRj ~ 12 (Table |6]l is required on 
the forest model and SNRm > 12 is achieved in a ~200-h obser- 
vation for the worst case of a G2 star at 5 pc (Fig. [19] top panel). 
SPICES can distinguish these three cases for any terrestrial plan- 
ets on a 1-AU orbit within 5 pc. If we consider 50% cloud cov- 
erage, surface effects are more difficult to distinguish and would 
require SNRi ~ 30 in the red part (Table |6]l. The exposure time 
exceeds the limit of 200 h for a target at 5 pc, which rather limits 
the sample to 4 pc (Fig. [19] bottom panel). Molecular oxygen 
and water absorptions as well as the "red edge" can still be mea- 
sured. On the contrary, 100% cloud coverage definitely prevents 
the identification of surfaces since it would require SNR, ~ 220. 
This performance that is out of reach of a small telescope like 
SPICES in a reasonable amount of time. In the favorable case of 
a Sun-like star at 2 pc (a Cen A is the sole known case), the in- 
strument achieves SNRm ~ 1 10 and allow to distinguish between 
cloudy planets totally covered with ocean and forest respectively. 

5. Potential targets 

In this section, we extend the results obtained above in order 
to determine the minimum planet radius accessible to SPICES 
for each planet category studied in the previous section. We also 
study the volume of the planet sample which can be character- 
ized (metallicity, cloud coverage, surface type) with SPICES. 
The models used are the gas giants with metallicity 1, the ice 
giants with metallicity 10, the cloudy forest Earths for the cloud 
coverage measurement, and the forest Earths with and 50% 
clouds for the surface type analysis. We consider that the planet 
flux increases as the square of its radius and that the atmosphere 
composition and structure remain the same. We set the following 
values of SNRri 30 for all Jupiters and Neptunes and 25, 12 and 
30 respectively for the Earths (Table |6). We assume a maximum 
exposure time of 200 h and three host stars of type G2, AO and 
MO. We carry out this study for the four models of Jupiter and 
Neptune at 0.8, 2, 5 and 10 AU and the Earth model at 1 AU. 
We recall that these separations are relevant for a G2 star and 
the corresponding separations for AO and MO stars are given in 
Table [5] We consider the following radius ranges for the planets; 



Rp < 2.5 Re for the Earths (iGrasset et al.ll2009l) . 
0.5 Rn <Rp < 1.5 Rn for the Neptunes. The lower limit is 
the radius for the maximum mass we consider for the Earths 
(10 Me). The upper value correspond s to a maximum mas s 
of 30 Me in the mass-radius relation of lFortnev et al.l(l2007l) . 
0.5 Rj <Rp < 1.1 Rj for the Jupiters. We derive the lower 
value from the upper mass limit we assume for the Neptunes. 
The upper limit corresponds to the maximum radius of a 4.5- 
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Fig. 20. Minimum radius of model planets for which character- 
ization is possible (see text) for different star types: G2 (top), 
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curves of the 0.8-AU Jupiter (top) and 0.6-AU Jupiter (bottom). 
Note that for lack of space, we do not indicate the labels for all 
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Gyr planet at 1 AU in iFortnev et aP (l2007h . We note that 
transit measurements find potentially larger planets (up to 
1 .4 Rj), but very close to the star. The inflated radius of these 
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Table 9. SPICES performance summary for exoplanet spectroscopic characterization. We strongly encourage the reader to refer to 
the rest of the paper to access the assumptions and the detail of the studies. 



Planet Probed property 



AO 



Stellar type 
G2 



Sections 



MO 



Jupiter 



lAU<a<4.2AU 0.25 AU<a<0.8 AU 



Rayleigh scattering 



Solar metallicity lx/3x 



5 pc <d< 17 pc 
0.5Rj<Rp<l.lR 



1.3pc<d<4pc 
0.5Rj<Rp<l.lR. 



1 AU < a < 4.2 AU 

5 pc <d< 17 pc 

0.6Rj<Rp<l.lRj 



0.25 AU<a<5 AU 

1.3pc<d<10pc 

0.5R|<R„<1.1R, 



0.25 AU < a < 3 AU 

1.3pc<d<7.5pc 

0.5R,<R„<1.1R, 



IIXni4:2l and[5] 



Neptune 



Rayleigh scattering 



Solar metallicity 10x/30x 



0.25 AU<a<0.8AU 

1.3pc<d<4pc 
0.5Rn<Rp<1.5Rn 

0.25 AU < a < 2 AU 

1.3pc<d<8.5pc 

0.5Rn<Rp<1.5Rn 



0.25 AU<a< 1.5 AU 

1.3 pc <d<3 pc 
0.5Rn<Rp<1.5Rn 



l2Ani43l and[5] 



Cloud coverage 0/50/100% 



Earth Forest coverage 0/50/100% 
clouds = 0% 



Forest coverage 0/50/100% 
clouds = 50% 



0.25 AU<a<l AU 

1.3 pc <d<5 pc 
1Re<Rp<2.5Re 

0.25 AU < a < 1 AU 

1.3 pc <d<5 pc 
1.1 Re<Rp<2.5Re 

0.25AU<a<l AU 

1.3pc<d<4pc 
1.1 Re<Rp<2.5Re 



0.25 AU < a < 0.6 AU 

1.3 pc <d<3 pc 
1.1 Re<Rp<2.5Re 

0.25 AU < a < 0.6 AU 

1.3 pc <d<3 pc 
1.2Re<Rp<2.5Re 

0.25 AU < a < 0.6 AU 

1.3pc<d<2.8pc 

1.4Re<Rp<2.5Re 



[2A2ll4:4l and[5] 



planets could be explained by add itional mechanisms to stel- 
lar radiation (iFortnev et al.ll2010l) . 

Figure |20] gives the minimum planet radius (in Re units) for 
which SPICES is able to characterize the metallicity of giants or 
the cloud and forest coverage of teiTestrial planets. The colors 
refer to the planet type (Jupiter, Neptune or Earth) and the line 
styles to the planet model (separation for the giants and cloud 
coverage for the Earths). 

For G2 stars (Fig. |20] top panel), SPICES reaches the low- 
est radius we consider (0.5 Rj ~5.5 Re) for the 0.8- and 2-AU 
Jupiters, for all stars within respectively 4 and 10 pc. About 
twenty G stars could be probed for Jupiters at 2 AU, while a few 
stars could be searched for Jupiters at 0.8 AU. For a Jupiter at 
5 AU, the upper limit in radius (1.1 Rj ~ 12 Re) is reached when 
the star is at 8.5 pc. The radius linearly decreases as the star dis- 
tance decreases down to 8 Re at 4 pc (flux proportional to the 
square ratio of the planet radius to the star distance). At shoiter 
distances, the deviations from the linearity are due to the speckle 
noise. When the star distance decreases, the angular separation 
of the star-planet system increases and the level of the speckle 
background increases, especially towards the edges of the dark 
hole (Fig.|6]l. As for Neptune-like planets, the minimum radius 
scales nearly linearly with distance from ~ 2 to 4 Re at 0.8 AU, 
and from 2 to 6 Re at 2 AU. SPICES can characterize cloudy 
and clear telluric planets at 1 AU around a few G2 stars up to 
5 pc (IWA limitation). Potentially, it is able to reach Earth-size 
planets for very close stars like a Cen A (~ 1.3 pc). 

For AO stars (middle panel), the only planet SPICES can 
study among the planetary atmosphere models we consider is 
the cloud-free Jupiter. We already noted the lack of models for 
separations smaller than ~4 AU (Sect. [33] ). The shoitest separa- 
tion available is 4.2 AU (Table |5]l, which is angularly resolved 
by SPICES when the star is closer than 21 pc. For this planet, 
the upper limit of the detectable radius roughly follows a linear 
law when the star distance increases from 9 to 17 pc. We note 



that a dozen of A stars are located within the upper limit of the 
range. For the closest stars, the detector read-out noise impacts 
the performance. 

If we focus on MO stars (bottom panel), we note that the 
cloud-free giant planets (0.24 AU, Table |5]l are inaccessible due 
to the IWA (Table[8l). For the Jupiters at 0.6 and 1 .5 AU (Table|5]l, 
we determine that SPICES allows the analysis of the smallest 
radius up to 3 and 4 pc respectively. The radius upper limit is 
achieved at ~7.5 and ~4 pc for separations of 1.5 and 3 AU 
respectively. Jupiters at 1.5 AU could be accessible for ~100 
M stars. For the 3-AU Jupiter case, the minimum radius that 
is detectable decreases as a linear function towards short dis- 
tances. The minimum radius (0.5-0.6 Rn) is feasible for 0.6-AU 
Neptunes (water clouds) at distances < 3 pc. SPICES can also 
access icy giants as far as 1.5 AU within ~3 pc (~10 M stars 
satisfy this constraint). Finally, telluric planets at 0.3 AU (the 
luminosity-scaled distance equivalent to 1 AU from a Sun-like 
star. Table |5]l are proven to be difficult to detect with SPICES 
and only reachable within 1.5 pc (Table |8]l. However, since M 
stars are of interest in terms of contrast, we extrapolate the lu- 
minosity of the model at 0.3 AU out to 0.6 AU, using a simple 
inverse square power law. This is obviously not rigorous, as the 
atmosphere characteristics would change, but it gives a rough 
estimation for one of the most challenging science cases of the 
mission. With such an assumption, we find SPICES can charac- 
terize telluric planets with radii as small as 1-2 Re within 3 pc. 

Table |9] summarizes SPICES' spectrophotometric perfor- 
mance for all planets and stellar types we considered in terms 
of star-planet separation (a), star distance (d), and planet ra- 
dius (Rp) ranges. As explained in Sect. 12.4.11 exoplanetary at- 
mosphere albedo can drastically change with the star-planet sep- 
aration. The minimum separation at which SPICES can detect 
a planet thus depends on the coronagraph IWA, the closest ob- 
servable star for each spectral type and the planetary atmosphere 
albedo. Because of the lack of atmosphere models for close-in 
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planets, we cannot derive a precise lower limit for the star- planet 
separation. However, we give a value that coiTesponds to the 
IWA limitation for the closest stars. For the star distance range, 
the lower limit is set to either the distance of the closest star 
of the considered spectral type, or the distance below which the 
planet is fainter than the speckle noise. The upper limit is either 
the maximum distance at which the planet separation is resolved 
by SPICES, or the distance beyond which the planet is fainter 
than the photon noise. 

6. Conclusions 

In this paper, we presented an analysis and an estimation of the 
performance of SPICES, a small coronagraphic space mission 
operating in the visible. From numerical simulations of the in- 
strument concept given in Boccaletti et al. (2012), we first de- 
termined that the contrast reached by the instrument meets the 
top-level requirements (-10^^ at 2 AID and ~10"'° at 4 AID). 
Then, we analyzed the impact of different sources of noise: de- 
tector read-out noise, exo-zodiacal intensity and photon noise 
considering several stellar types. We confirmed previous analy- 
ses that exodisks might be a major limitation for the characteri- 
zation of faint planets, like Neptunes and super-Earths, because 
it biases the planet photometry (for disk intensity > 1 zodi) and 
adds photon noise (>a few 10 zodis). Exodisk intensity calibra- 
tion and target selection will be necessary to reduce this limi- 
tation. Using planet models calculated for a Sun-like star and 
assuming flux conservation, we found that the brightest stars 
(types AF) would not be the most favorable targets for SPICES. 
For instance, only planets with clear atmospheres will be de- 
tectable around AO stars within ~20 pc. On the contrary, stars of 
types GKM will allow the detection of a large variety of plan- 
ets as far as ~ 12 pc for a G2 star and ~7.5 pc for a MO star: 
cloud-free, water-cloud and ammonia-cloud Jupiters, cloud-free 
and water-cloud Neptunes and telluric planets. After this gen- 
eral study of contrast performance, we focused on spectrometiic 
abilities on planets around a solar-type star. We defined a crite- 
rion on the SNR of the measured spectra to determine if SPICES 
or an analog space mission could probe several planetary prop- 
erties by disentangling their spectra (clouds and metallicity for 
giant planets, cloud and surface coverage for telluric planets). 
The results are very encouraging since SPICES could character- 
ize the metallicity of Jupiters and Neptunes up to 10 pc and 6 pc 
respectively for a G2 star. The cloud and surface coverage of a 
super-Earth orbiting a Sun analog would also be characterized 
for systems closer than 5 pc. For AO stars, the instrument could 
study the metallicity of cloud-free Jupiters closer than 17 pc. For 
MO stars, cloudy Jupiters and Neptunes would be accessible at 
distances <7.5 pc and <3 pc respectively, and super-Earths at 
distances <3 pc. These results would give a potential target list 
of ~300 stars. We also estimated the minimum radius for the 
planets SPICES could detect. 

We emphasize that our study is also useful for other space 
coronagraph studies cuiTently under development (Sect. [!]• 
Indeed, as far as the authors of this paper are aware of, no study 
has investigated in detail the spectrophotometiic abilities of such 
missions to retrieve physical parameters of realistic planet spec- 
tra. An important point we noted during our study is the need 
for grids of planetary spectra, as it has been done for NIR 
planet-finders on 8-10 m telescopes (e.g ., Chabrier et al. 2000; 
iBaraffe et all l2003t iBurrows et all l2006t iFortnev et al.l |2008|). 
For self-luminous planets, the problem is a bit simplified with 
respect to mature planets, because the emergent spectrum is in- 
dependent of the separation from the star and the stellar irradi- 



ation for separations > 1 AU. Reflected spectral models are usu- 
ally limited to solar- type stars and are derived for a few points of 
the planet parameter space (mass, separation, metallicity, age). 
To thoroughly estimate instrument performance, we will need 
spectra for large ranges of planet types (gas and ice giants, and 
super-Eaiths), separations, stellar types, etc. 

In this paper, we also consider areas for future work that will 
refine our results. We will include Fresnel propagation in our in- 
strument model to study the impact of out-of pupil aberrations 
on the performance. Phase and amplitude aberrations will evolve 
with the wavelength. Their calibration will be as good as they 
cuiTently are, because our focal plane wavefront sensor can esti- 
mate both phase and amplitude abeiTations in each channel of the 
IFS independently. The DM correction may be slightly degraded 
and an optimization of the design may be required to reach a 
contrast of 10"'". We also plan to study the polarimetric perfor- 
mance of SPICES, as we expect that the combination of flux 
and polarization measurements can remove degeneracies that 
arise when retrieving planet properties from flux measurements 
alone. Circumstellar disks are another science case to study 
(dust distribution, rings, planet gaps). Finally, another point of 
great interest would be to apply the methods we used here to a 
larger telescope lik e the Terrestrial Planet Finder Coronagraph 
dLevine et alj|2009l) . Indeed, one of the main results of our work 
is that the optimal targets of SPICES would be at separations 
of ~0.8-2 AU, because of the drastic decrease of the reflected 
flux. For these cases, the performance is limited by the angular 
resolution of the telescope and not the photon noise. The choice 
of the telescope diameter was anterior to this study and resulted 
from a trade-off between the science objectives and the budget 
allocated to a medium-class mission. We now consider to sub- 
mit a larger (2.5-3 m) instrument as a large-class mission when 
a call for proposals will be issued. Besides these considerations, 
such a telescope will access a larger volume of target stars and 
will be less sensitive to zodiacal and exo-zodiacal contributions 
for planet spectra measurement (Traub & Oppenheimer2010) . 
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